function Htot=Mod_Filtros_Mix(H1,H2,transicional)

   

    [zb pb kb]=zpkdata(H1);
    [zc pc kc]=zpkdata(H2);
    pb=cell2mat(pb);
    pc=cell2mat(pc);
    zb=cell2mat(zb);
    zc=cell2mat(zc);
    for k=1:length(pb)
        if abs(imag(pb(k)))<abs(pb(k))/1000
            pb(k)=real(pb(k));
        end
        if abs(imag(pc(k)))<abs(pc(k))/1000
            pc(k)=real(pc(k));
        end
    end
    for k=1:length(pb)
        pb(k)=-abs(real(pb(k)))+j*imag(pb(k));
        pc(k)=-abs(real(pc(k)))+j*imag(pc(k));    
    end
    for k=1:length(pb)
        for m=k:length(pb)
            if (imag(pb(m))/abs(pb(m)))>(imag(pb(k))/abs(pb(k)))
                aux=pb(m);
                pb(m)=pb(k);
                pb(k)=aux;
            end
            if (imag(pc(m))/abs(pc(m)))>(imag(pc(k))/abs(pc(k)))
                aux=pc(m);
                pc(m)=pc(k);
                pc(k)=aux;
            end


        end
    end

    for k=1:length(pb)
        ptot(k)=pb(k)*(1-transicional)+pc(k)*transicional;
    end
    
    if isempty(zb)
       ztot=[];
    else
       for k=1:length(zb)
            for m=k:length(zb)
                if imag(zb(k))>imag(zb(m))
                    aux=zb(m);
                    zb(m)=zb(k);
                    zb(k)=aux;
                end
                if imag(zc(k))>imag(zc(m))
                    aux=zc(m);
                    zc(m)=zc(k);
                    zc(k)=aux;
                end


            end
       end
       for k=1:length(zb)
       ztot(k)=zb(k)*(1-transicional^1.1)+zc(k)*transicional^1;
    end
        
    end
    

    Htot=tf(poly(ztot),poly(ptot));
    mag=bode(Htot,1e-5);
    Htot=Htot/mag;


end